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FINAL  REPORT 


DYNAMIC  STRUCTURE  LES  MODELS 

AFOSR  GRANT  NUMBER:  F49620-02-1-0348 

Christopher  J.  Rutland 
Engine  Research  Center 
Department  of  Mechanical  Engineering 
University  of  Wisconsin  -  Madison 

Abstract 

A  new  dynamic  structure  (DS)  approach  for  large  eddy  simulation  modeling  has  been 
developed.  The  approach  uses  a  tensor  coefficient  for  the  structure  of  the  sub-grid  stresses 
and  a  sub-grid  kinetic  energy  for  the  magnitude  of  the  stresses.  The  tensor  coefficient  is 
obtained  directly  from  using  the  common  ‘dynamic  approach’  using  the  Germano  identity. 
The  sub-grid  kinetic  energy  is  obtained  from  a  transport  equation.  The  objectives  of  the 
program  are:  (i)  to  develop  the  DS  approach  for  scalar  transport  with  potential  application 
to  heat  transfer  problems  and  (ii)  develop  and  test  the  DS  formulation  for  rotating 
turbulence.  The  approach  offers  several  significant  advantages  over  existing  approaches 
including:  (i)  no  turbulent  viscosity  is  used,  (ii)  an  energy  budget  is  enforced,  (iii)  very 
good  prediction  of  the  actual  sub-grid  stresses  is  obtained,  and  (iv)  the  model  rests  on  a 
sound  mathematical  foundation  that  observes  solvability  constraints.  New  subgrid  scale 
models  are  developed  for  subgrid  scalar  flux,  subgrid  scalar  dissipation  and  subgrid  scale 
energy  dissipation.  Models  are  evaluated  a-priori  using  available  DNS  data  for  decaying 
isotropic  turbulence,  incompressible  channel  and  Couette  flows,  and  non-reacting  mixing 
layer.  A  posteriori  tests  include  LES  of  decaying  isotropic  turbulence  and  LES  of 
incompressible  mixing  layer. 

Background 

The  next  generation  of  turbulent  flow  simulations  will  undoubtedly  use  large  eddy 
simulation  (LES)  modeling  approaches  (Moin  1997,  Ferziger  1983).  However,  the 
adoption  of  LES  for  more  general  use  has  been  slow,  in  large  part  due  to  the  very  large, 
nearly  DNS,  grids  that  are  required.  One  could  attribute  the  large  grid  size  requirements  to 
limitations  of  the  current  modeling  approaches.  Recently,  a  new  approach,  called  dynamic 
structure  modeling  (Pomraning  2000,  Pomraning  and  Rutland  2002),  has  been  developed 
that  offers  several  important  advantages  over  existing  LES  approaches. 

In  LES,  spatial  filtering  is  used  to  remove  scales  that  must  be  resolved  in  the  simulation. 
When  applied  to  incompressible,  non-reacting  turbulent  flow,  the  basic  term  that  requires 
modeling  for  is  the  sub-grid  stress  tensor: 

(1) 

where  the  over-bar  indicates  LES  spatial  filtering.  The  dynamic  Smagorinsky  approach  is 
the  most  commonly  used  for  current  applications.  This  model  assumes  that  the  stress 
tensor,  Ty,  can  be  linearly  related  to  the  resolved  scale  strain  rate  tensor: 

20050309  085 


1 


(2) 


where  StJ  is  the  resolved  scale  strain  rate  tensor  and  vr  is  a  ‘turbulent’  viscosity 

commonly  modeled  as  vT  =  Cs A2  |sj  with  A  a  representative  length  scale  and  Cs  a 

modeling  coefficient.  (Note,  for  clarity  a  trace  term,  commonly  absorbed  into  the 
incompressible  pressure,  is  not  shown.)  This  is  a  turbulent  viscosity  based  model,  and  any 
model  based  on  a  viscosity  approach  will  require  a  positive  viscosity  coefficient  for 
numerical  stability.  This  means  that  the  sub-grid  model  is  only  a  dissipative  term  and 
removes  kinetic  energy  from  the  flow.  However,  Ty  is  not  a  purely  dissipative  term  -  it  is 

a  local  term,  because  LES  filtering  is  a  local  operation.  Moreover,  the  application  of  eddy 
viscosity  models  in  LES  violates  one  of  the  basic  assumptions  of  gradient  modeling  -  that 
the  length  scale  of  the  modeled  phenomena  is  significantly  smaller  than  the  distance  over 
which  the  mean  gradient  changes  appreciably,  as  pointed  out  by  Corrsin  (Corrsin  1974). 


The  dynamic  approach  uses  another  ‘test’  level  filtering  operation  (indicated  by  a  peaked 
over-bar)  and  a  new  tensor  arising  from  the  Germano  identity  (Germano  1991): 

Ly  =uiUj-uiuj  (3) 

This  tensor  is  known  from  the  resolved  scale  and  can  be  used  to  formulate  the  following 
equation  for  the  model  coefficient  in  the  dynamic  Smagorinsky  model: 

L,-\sju=afi,-P&  (4) 


where  a,.  =  -2A 


StJ  and  /3tj  =  -2A2|s|S(/  are  introduced  for  notational  simplicity.  This 


is  an  extra  equation  that  can  be  used  to  solve  for  either  the  viscosity  directly  or  the  model 
coefficient  in  the  viscosity  of  Eq.  (2).  However,  Eq.  (4)  is  a  second  rank  tensor  equation 
(e.g.  six  equations  by  symmetry)  for  a  single  unknown  coefficient,  Cs .  This  highly  over 

specified  situation  is  commonly  dealt  with  by  using  a  minimization  procedure  to  find  the 
coefficient.  We  note  that  this  over  specification  could  be  interpreted  as  a  fairly 
fundamental  problem  with  the  whole  dynamic  Smagorinsky  approach. 


There  are  two  other  approaches  in  LES  modeling  due  to  Ghosal  et  al.  (1995)  and  Menon  et 
al.  (1995).  Both  of  these  use  a  sub-grid  kinetic  energy  and  this  represents  an  advance  over 
the  zero  equation  Smagorinsky  approach.  However,  both  approaches  use  the  kinetic 
energy  only  as  an  additional  scaling  term  to  form  a  turbulent  viscosity.  Thus,  they  still 
follow  Eq.  (2)  and  suffer  from  the  same  inherent  limitations. 


The  Dynamic  Structure  Approach 


To  address  some  of  the  mathematical  and  conceptual  difficulties  in  the  common  LES 
models,  a  new  modeling  approach  was  developed.  This  new  approach  does  not  use  a 
turbulent  viscosity  concept,  but  instead  attempts  to  model  the  sub-grid  stresses  directly. 
This  is  accomplished  by  using  the  following  model  where  k  is  the  subgrid  kinetic  energy: 


h 


=  cyk 


(5) 
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and  Cy  is  a  tensor  coefficient  that  is  obtained  directly  as  part  of  the  solution  using  the 
dynamic  procedure.  This  is  accomplished  by  defining  a  ‘test’  level  model: 

Ty=CyK  K  =  -«'«’)  =  k+^L„  (6) 

using  the  test  level  kinetic  energy,  K ,  can  be  written  as  a  function  of  the  sub-grid  kinetic 
energy  and  the  trace  of  the  Ly  tensor.  Then  the  Germano  identity  is  used  to  formulate  an 
equation  for  the  tensor  coefficient,  ctJ : 

Ly=KCy~Uy  (7) 

Note  that  this  is  a  well-poised  problem  in  that  the  number  of  coefficient  unknowns 
matches  the  number  of  equations.  Fundamentally  this  is  an  integral  equation  for  the 
coefficient  tensor  (e.g.  the  over-bar  indicates  an  integral  filter)  and  we  have  tested  it  in  the 
integral  form  (Pomraning  2002).  However,  a  very  effective  algebraic  form  of  the  model  is 
formulated  by  assuming  that  the  coefficient  tensor  can  be  brought  outside  the  integral. 
Then,  the  dynamic  structure  model  becomes: 

T«=(V4i)2t  (*) 

Here  the  sub-grid  kinetic  energy  is  obtained  from  a  transport  equation.  Thus,  the  dynamic 
structure  model  is  a  one-equation,  non-turbulent  viscosity  model.  The  tensor  structure  of 
the  sub-grid  stress  tensor,  xtj ,  is  found  directly  from  the  dynamic  approach  through  the 

tensor  coefficient.  This  is  very  advantageous  because  the  coefficient  is  just  the  normalized 
Ltj  tensor,  which  is  directly  related  to  the  structure  of  the  sub-grid  stresses  (see  Eq.  (3)). 1 

Then,  the  magnitude  of  is  obtained  from  the  sub-grid  kinetic  energy,  k,  which  is 

simply  twice  the  trace  of  ry  itself. 

Advantages  of  the  New  Approach 

1.  Non-viscosity  based  -  By  not  using  a  turbulent  viscosity  approach,  the  dynamic 
structure  model  is  not  purely  dissipative.  This  allows  local  solutions  that  respond  to  the 
local  flow  conditions  -  conditions  that  are  not  always  dissipative.  This  suggests  that  the 
model  has  the  potential  to  give  the  correct  behavior  in  flows  with  streamline  curvature 
or  system  rotation  that  often  cause  traditional  approaches  difficulty.  There  is  no 
turbulent  viscosity  in  the  new  model  and  molecular  viscosity  plays  its  correct  physical 
role  of  dissipating  sub-grid  energy. 

2.  An  energy  budget  -  The  dynamic  structure  model  has  a  sub-grid  transport  equation  for 
k .  Hence  it  is  a  one-equation  model  rather  than  a  zero  equation  model  like  the  dynamic 
Smagorinsky  model.  This  energy  budget  keeps  the  model  stable  without  a  turbulent 
viscosity.  Hence,  it  is  very  robust  and  works  over  a  wide  range  of  grid  resolutions 
including  practical  situations  in  which  there  is  significant  energy  in  the  sub-grid  scales. 


1  Note  that  Meneveau  group  (Meneveau  and  Katz  2000,  Liu  et  al.  1998)  once  suggested  a 
similar  idea  but  did  not  find  the  appropriate  scaling  for  the  magnitude  of  the  sub-grid 
stress  tensor. 
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Dynamic  Smagorinsky  models  require  dense  grids  with  little  energy  in  the  sub-grid 
scales  so  that  the  sub-grid  model  is  not  required  to  do  much. 

3.  Good  prediction  of  the  sub-grid  stresses  -  The  component  structure,  spatial  distribution, 
and  magnitude  of  the  sub-grid  stresses,  zv ,  are  very  well  predicted  by  the  new  model. 
This  is  shown  in  Fig.  1  that  compares  the  new  model  to  the  standard  dynamic 
Smagorinsky  model.  Clearly  the  new  model  does  a  much  better  job  of  predicting  the 
actual  sub-grid  stresses.  This  is  due  to  the  tensor  coefficient  that  comes  directly  from 
using  the  stress  tensor  structure  in  the  dynamic  approach.  This  makes  the  new  model  a 
very  strong  candidate  for  flows  with  more  complex  physics  such  as  heat  transfer 
because  the  sub-grid  terms  can  be  used  to  build  other  sub-models. 

4.  Solvability  -  The  mathematical  formulation  of  the  new  model  is  much  sounder  than  the 
common  dynamic  Smagorinsky  model.  The  underlying  formulation  of  the  dynamic 
approach  is  a  Fredholm  integral  equation  of  the  second  kindjhat  has  specific  solvability 
constraints.  Other  LES  formulations  that  relate  xtj  to  Sy  can  be  shown  to  violate 
solvability  criteria  for  almost  any  type  of  flow.  However,  the  dynamics  structure  model 
relates  Ti}  to  the  sub-grid  kinetic  energy  and  this  results  in  integral  equations  that  easily 
satisfy  solvability.  This  sound  mathematical  formulation  adds  robustness  to  the  model 
so  that  it  works  under  a  wide  variety  of  flows  and  grid  resolutions.  Additionally,  it  also 
appears  to  guarantee  other  nice  properties  such  as  realizability  (e.g.  non-negative 
diagonal  stress  components)  .  Finally,  we  note  that  the  model  is  Galilean  invariant 
(Pomraning  2002). 


Figure  1:  Results  from  isotropic  decaying  turbulence;  contours  of  the  r12  sub-grid  shear 
stress:  (a)  DNS  result,  (b)  dynamic  structure  model,  (c)  dynamic  Smagorinsky  model 


The  dynamic  structure  model  (Eq.  (8))  has  been  tested  in  a  high  order  LES  code  for 
isotropic  turbulence  and  compared  very  well  with  Comte-Bellot  results  (Pomraning  2002). 
Note,  unlike  many  Smagorinsky  based  models,  the  sub-grid  kinetic  energy  accounted  for  a 
significant  portion  of  the  total  kinetic  energy,  often  1 0%  or  more. 

The  model  has  also  been  applied  to  other  problems  in  a  low-order  ‘engineering’  code. 
While  there  are  issues  of  numerical  accuracy  in  these  codes,  this  is  somewhat  mitigated  by 
the  higher  Reynolds  numbers  and  courser  grid  resolutions  that  require  larger  magnitudes 
of  the  model  terms.  Some  results  are  shown  for  the  Eaton  et  al.  (1986)  backwards-facing 
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step  in  Figure  2.  These  results  required  time  averaging  of  the  LES  results  to  achieve  the 
appropriate  comparison  with  the  time  averaged  experimental  results.  This  figure  also 
shows  the  results  of  using  no  model  to  demonstrate  that  the  numerical  dissipation  is  not 
providing  the  correct  sub-grid  model. 


-1.0  -0.5  0.0  0.5  1.0  1.5  -1.0  -0.5  0.0  0.5  1.0  1.5  -1.0  -0.5  0.0  0.5  1.0  1.5 


u/U  u/U  u/U 

Figure  2:  Normalized  streamwise  velocity  profiles  at  three  locations  downstream  of  a 
backwards  facing  step;  comparisons  between  experimental  data  of  Eaton  et  al.  [9],  a 
k-e  case,  the  dynamic  structure  model,  and  a  no-model  (e.g.  ‘laminar’)  case. 


LES  Equations  for  Scalar  Transport 


The  LES  transport  equation  for  a  scalar  <f>  is  obtained  by  applying  the  spatial  filtering 
operation  to  the  original  transport  equation  to  yield: 

dl  t  duj*_  d  La^l  3tm 

dt  dxj  dxj  dxj  dxj 


Here,  ^denotes  the  spatially  averaged  scalar  D  is  the  diffusion  coefficient,  and 
tit  =  Vfi-V~4>  is  the  subgrid-scale  (SGS)  flux  term  which  has  to  be  modeled.  The  most 
popular  models  are  gradient-type  models: 


Various  techniques  for  determination  of  the  turbulent  viscosity  Hj  have  been  proposed 
(Kim  et.al  1999;  Kimmel  and  Domaradzki  2000).  The  dynamic  viscosity  approach  uses 
another  filtering  operation  referred  to  as  “test”  filtering  and  denoted  by  1j>  (Germano  et.al. 
1991).  Similar  to  models  for  SGS  Reynolds  stress,  the  eddy  viscosity  models  are  hardly 
applicable  for  turbulence  modeling,  as  outlined  by  Corrsin  (1974).  Another  approach  is  to 
use  self-similarity  models  of  the  type 

TjAaC\_uj~uj^\  (11) 
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where  C  is  a  scaling  constant,  or  mixed  models  that  combine  (10)  and  (1 1)  (Speziale  et.al. 
(1988)) 


Dynamic  Structure  Models  for  SGS  Scalar  Flux  and  Scalar  Dissipation 


In  the  spirit  of  the  model  (8)  developed  by  Pomraning  and  Rutland  (2002),  we  propose  the 
following  model  for  the  subgrid  scalar  flux: 


(12) 


Here  =  is  the  SGS  variance  of  <j>  on  the  base  (grid) 


level,  0=<#<*-^is  the  SGS 


variance  of  <j> on  the  test  level,  and  Ltj)  =  u<j>  -  ut<j> is  the  modified  Leonard-type  term  for 

scalar  equation.  Note  that  the  SGS  variance  9  cannot  be  obtained  from  the  resolved 
quantities  and  thus  has  to  be  modeled  (Cook  and  Riley,  1994;  Pierce  and  Moin;  1998)  or 
obtained  from  a  separate  transport  equation  (Jimenez  et.al,  2001).  In  keeping  with  the  DS 
approach,  we  use  a  transport  equation  for  6  that  can  be  obtained  from  (9): 


80  8Uj0 

dt  dxj 


_a_ 

dx. 


dx. 


-  dr. . 


(13) 


where 

and  C,  =-^-|Vu  -tfui 

are  the  SGS  scalar  dissipation  term  and  triple  correlation  term,  respectively.  The  triple 
correlation  term  is  usually  either  omitted  or  modeled  using  gradient-type  models. 
Alternatively,  it  is  possible  to  approximate  the  sum  of  last  two  terms  in  the  transport 
equation  for  9  to  the  third  order  of  computational  mesh  spacing  using  the  series  given  by 
Bedford  and  Yeo  (1993).  The  resulting  equation  is: 


X,=2D 


d<t>  d<f>  d<j>  d<j> 
dXj  dXj  dXj  dxj 


SO  dufi 

dt  dx. 


dx,. 


D 


de_ 

dx 


7J 


iA 


(14) 


where  the  subgrid  flux  is  modeled  using  (12)  and  the  SGS  scalar  dissipation  is  modeled 
using  its  own  modified  Leonard  type  term: 


xt 


,0,  .e_  did£ 

©  1  0  dxj  dXj 


dt_8l_ 

dXj  dxj 


(13) 


A  priori  Test 

The  scalar  flux  and  scalar  dissipation  models  were  tested  using  a  DNS  data  for  non¬ 
reacting  incompressible  spatially  developing  mixing  layer  from  Mason  and  Rutland 
(2000).  The  Reynolds  number  based  on  the  inlet  vorticity  thickness  was  200.  For  the 
purpose  of  the  a  priori  testing,  the  modeled  quantities  were  computed  directly  from  the 
DNS  data  using  filtered  DNS  results  and  the  LES  models.  In  order  to  test  robustness  of 
the  models,  LES  filters  of  various  dimensions  were  applied  to  the  DNS  data.  The  DS 
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models  for  the  SGS  scalar  flux  and  SGS  scalar  dissipation  were  compared  to  several 
models  in  the  literature. 

Figure  3a  represents  the  results  of  comparison  of  the  DS  model  (12)  and  a  viscosity 
gradient  model  (10).  The  DS  model  has  a  smaller  data  scatter  indicating  that  it  performs 
better.  These  comparisons  were  made  over  a  range  of  LES  filter  sizes  and  shapes  and  the 
DS  model  performed  better  in  all  cases  (Chumakov  and  Rutland,  2003) 


Figure  3:  (a)  Scatter  plot  of  SGS  scalar  flux  term  computed  from  DNS  and  using  DS 
model  (blue  circles)  and  viscosity  models  (red  triangles)  (b)  PDF  of  relative  error  for 
models  of  the  SGS  dissipation  term:  DS  model  (blue  solid)  and  Jimenez  et  al.  (2001)  (red 
dashed).  Different  curves  represent  different  LES  filters  shapes  and  sizes. 


Figure  3b  depicts  the  results  of  a  similar  a  priori  test  for  two  SGS  dissipation  models.  The 
DS  model  (12)  is  compared  to  the  one  from  Jimenez  et.al  (2001) 


“  dXj  dXj  k 


where  s0  =  v(du,  /dxj  du,  Idxj)  is  the  SGS  kinetic  energy  dissipation.  Note  that  this  model  is 
similar  to  a  RANS  type  models  in  that  it  uses  a  turbulence  time  scale  for  the  scalar 
dissipation.  Experience  has  shown  that  this  time  scale  is  often  not  appropriate  for  scalar 
quantities,  for  example  in  turbulent  regions  where  the  scalar  quantity  has  become  uniform 
from  mixing.  Figure  lb  shows  a  PDF  plot  of  the  a  priori  error  for  the  scalar  dissipation 
models.  The  DS  model  performs  better  with  errors  centered  at  zero  and  smaller  spread. 


Subgrid  Scale  Energy  Dissipation 


The  SGS  kinetic  energy  is  obtained  from  a  transport  equation,  which  requires  a  model  for 
SGS  energy  dissipation 


€SGS 


duf  dut  diii  du^ 
dx.  dx.  dx.  dx. 


Together  with  the  energy  transfer  term  n  =  ry  Sy ,  the  SGS  dissipation  plays  a  crucial  role 

in  the  energy  transfer  process  between  the  turbulence  scales.  While  II  is  responsible  for 
the  energy  transfer  between  the  resolved  and  unresolved  scales,  s  removes  the  energy 
from  the  unresolved  scales.  A  model  has  been  tested  for  e  which  has  the 
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form  fsoj  »  v-F-Lc  where  Le  is  the  Leonard  type  term  associated  with  dissipation  (see  Eq. 
(14)).  The  scaling  factor  F  was  investigated  from  a  priori  tests.  Three  DNS  datasets  of 
decaying  isotropic  with  ReA  « 104.5,  92.6,  and  53.7  were  used.  The  proposed  form  of  the 

scaling  factor  F  is:  f = C,  [k/K] ,  where  2 K  =  L„ =u^-u^  is  the  trace  of  the  Leonard  term, 
and  Cc  and  y  are  the  constants  to  be  determined  from  the  a  priori  study.  Probability 
density  functions  (PDFs)  for  Cc  show  a  good  collapse  for  y  =  1/2 .  An  additional 
dependence  on  the  base  filter  to  the  3/2  power  was  also  seen.  Thus,  a  length  scale  ratio 
was  added  to  F  so  that  the  form  of  the  model  is: 

£-*K'c[y]  S 

where,  C  is  a  scaling  constant  and  /  is  a  base  filter  level  length  scale,  both  to  be 
determined  in  future  work. 

Limiting  Behavior  of  Dynamic  Structure  Models 

Dynamic  structure  models  use  ratio  of  variances  as  a  scaling  factor.  We  examine  the 
limiting  behavior  of  the  scaling  factors  <9/0  and  2 k/Lmm  for  the  case  when  the 

denominator  of  the  scaling  factor  goes  to  zero.  This  can  occur  when  either  the  grid  is 
asymptotically  refined  (A  ->0),  or  the  flow  field  becomes  quiescent,  for  example  in  the 
far-field  of  a  jet.  For  simplicity,  only  <9/0  is  presented.  Using  a  Leonard  expansion 
(Bedford  and  Yeo  (1993),  Leonard  (1974))  one  can  show  that  for  both  cases  the  scaling 
factor  stays  bounded  by  minimum  and  maximum  of  ratios  of  directional  second  moments 
of  the  base  and  test  filters,  i.e., 

ot  0  cc 

0<  min—  <  —  <  max—  <qo,  where  at  =  I  x~G(x)  dxi ,  a,  =  I  x;  G(x)  dxi . 

a  ©  a  i  -i 

This  statement  is  illustrated  for  a  one-dimensional  scalar  function  <j> .  Figure  4  presents  a 
scatter  plot  of  the  mean  base  variance  6  vs.  the  mean  scaling  factor  <9/0  together  with 
the  line  y  =  ata .  The  plot  was  obtained  by  successively  filtering  a  random  signal  to 
obtain  increasing  smoother  signals  (and  hence  smaller  variances).  This  shows  that  the 

mean  value  of  the  scaling  factor  goes  to  the  ratio  of  a! a  as  the  scalar  field  becomes 
quiescent  as  indicated  by  the  diminishing  SGS  variance. 


8u,  du,  du,  dll. 


dXj  dXj  dxj  dXj 


(14) 
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Figure  4.  Scatter  plot  of  the 
mean  base  SGS  variance  vs.  the 
mean  scaling  factor,#/©.  The 

line  represents  the  ratio  at  a 
and  indicates  the  limiting 
behavior. 


A-posteriori  Test  -  non-reacting  incompressible  mixing  layer 

In  order  to  evaluate  both  DS  models  an  LES  simulation  of  a  non-reacting  incompressible 
spatially  developing  mixing  layer  was  performed.  All  parameters  in  the  LES  simulation 
were  as  close  to  DNS  parameters  as  possible  but  the  LES  used  V*  as  many  grid  points  in 
each  coordinate  direction.  The  first-  and  second-order  statistics  were  computed  at  various 
locations  downstream  and  then  compared. 


Figure  5:  (a)  Rates  of  growth  for  mixing  layer  based  on  1%  scalar  difference:  DNS  and 
LES;  (b)  Comparison  of  the  scalar  variance  for  DNS  and  LES. 


Figure  5a  shows  the  comparison  of  first-order  statistics  from  the  DNS  and  LES  runs.  Both 
DNS  and  LES  mixing  layer  exhibit  linear  rate  of  growth  in  the  domain,  as  generally  found 
in  literature  (Dimotakis,  1991).  The  growth  rate  was  determined  using  the  1%  passive 
scalar  thickness  and  a  velocity  ratio  of  1/3.  The  growth  rate  coefficient  for  the  LES  was 
found  to  be  0.455.  This  is  similar  to  the  DNS  value  of  0.422  and  close  to  the 
experimentally  accepted  range  of  0.25  to  0.45.  Figure  4b  shows  the  comparison  between 
the  LES  and  DNS  mean  scalar  profiles  near  the  end  of  the  computational  domain. 
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The  variance  of  the  passive  scalar  from  the  DNS  and  LES  simulations  is  presented  in 
Figure  5b.  There  is  good  agreement  in  magnitude  and  spanwise  distribution.  Note  that 
there  is  a  large  contribution  from  the  subgrid  variance  9  to  the  total  variance  of  the  scalar 
-  a  characteristic  common  of  DS  type  models. 

A-posteriori  Test  -  decaying  isotropic  turbulence 

LES  of  decaying  isotropic  turbulence  is  a  good  test  for  both  the  SGS  stress,  ,  model  and 
the  SGS  energy  dissipation,  eSGS,  model.  In  general,  decaying  turbulence  exhibits  two 

different  types  of  behavior:  an  inertia-dominated  earlier  stage,  and  viscosity-dominated 
(final)  stage  (Pope  2000).  The  energy  decay  rate  for  the  final  period  can  be  obtained 
analytically  from  the  Karman-Howarth  equation,  which  yields  E  ~  t~sn .  For  the  inertia- 
dominated  period,  only  experimental  measurements  are  known,  but  the  total  kinetic  energy 
clearly  exhibits  the  power-law  decay  suggested  by  Comte-Bellot  and  Corrsin  (1966): 
E~r" .  The  values  of  the  decay  exponent  n  between  1.15  and  1.45  are  reported  in  the 
literature  but  Mohammed  and  LaRue  (1990)  suggested  that  nearly  all  of  the  data  are 
consistent  with  n  =  1 .3 .  Thus  we  match  the  energy  decay  rate  given  by  LES  against  the 
power  law  with  n  =  1 .3  for  the  inertia-dominated  period,  and  n  =  2.5  for  the  final  period. 


Figure  6.  Decay  of  total  kinetic 
energy  in  time.  Slopes  for  n=1.3 
and  n=2.5  are  plotted.  (Initial 
Taylor  scale  Reynolds  number  = 
200,  grid  size  323.) 


1  10  100  1000 
The  dynamic  structure  model  and  another  ^-equation  LES  model  are  compared.  The  other 
model  is  a  viscosity  based  model  and  uses  tv  w  -0.05-JkASijy  e  »1.0&3/2/A  (“Localization 

models”,  LM,  Ghosal  et  al.  1995,  Kim  and  Menon  1995,  Menon  et.  al,  1996).  Figure  6 
shows  the  log-log  plot  of  the  decay  of  total  kinetic  energy  in  time  for  both  models,  along 
with  the  slopes  for  n  =  1.3  and  n  =  2.5 .  In  the  LM  run,  the  transition  to  the  final  phase  of 
the  flow  is  more  spread  out  in  time  -  approximately  between  t=20...60,  as  opposed  to 
between  t=30...50  for  the  DS  run.  The  main  difference  between  the  two  models  is  the  ratio 
of  energy  in  the  subgrid  scales  to  the  total  kinetic  energy.  In  the  latter  stages  of  decay  this 
ratio  should  diminish.  This  was  seen  only  in  the  DS  model.  The  LM  model  showed  a 
nearly  constant  ratio  of  ~0.1  throughout  the  simulation  (Chumakov  and  Rutland  2005). 

Rotating  Flows 

The  dynamic  structure  approach  is  being  extended  to  simulate  rotating  flows.  Since  the  DS 
approach  maintains  an  energy  budget  between  resolved  and  subgrid  scales,  it  is  a  good 
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candidate  for  modeling  rotating  flows  in  which  there  can  be  strong  energy  transfer  to 
larger  scales.  The  subgrid  k  equation  in  the  DS  approach  does  not  require  modification 
since  the  Coriolis  term  only  appears  in  the  momentum  equation.  Thus,  efforts  are  focused 
on  studying  the  subgrid  stresses  in  rotating  systems. 


In  a  rotating  frame  of  reference,  the  subgrid  stresses  are  modified.  Models  for  the  subgrid 
stresses  should  also  be  modified  to  maintain  a  material  frame  indifference,  MFI,  (see 
Speziale,  1985,  and  Shimomura,  1999).  Two  models  which  observe  MFI  are: 


(15) 


where: 


^ ij-UiUj  UiUj  +  Cc 


(16) 


du,  duj 
dxk  dxk 


(17) 


Both  of  these  models  use  the  dynamic  structure  methodology  of  using  a  subgrid  kinetic 
energy  to  scale  a  normalized  tensor.  In  equation  (16)  the  first  terms  on  the  right-hand-side 
are  the  modified  Leonard  terms  that  appear  in  the  original  DS  model.  The  additional  terms 
in  brackets  represent  a  cross  correlation  term  is  required  to  maintain  MFI.  The  cross 
correlation  terms  are  modeled  using  an  additional  filtering  level  and  a  coefficient,  Cc ,  of 
0(1)  (set  tol.5  below).  In  the  second  model  (17),  the  nonlinear  term  Gi}  is  derived  by 
Taylor  expansion.  It  is  tested  because  it  maintains  MFI. 


Two  correlation  coefficients,  defined  in  Figure  7,  are  used  for  a  priori  testing:  The 
first,  Pj,  is  the  correlation  coefficient  (e.g.  the  scatter),  and  the  second,  p2,  is  the 
regression  coefficient  (e.g.  the  slope  of  the  scatter  line).  Each  correlation  goes  to  unity 
when  a  and  b  are  perfectly  correlated. 


_  <ab>-<axb> 

Px  <*/(<  a2  >  -  <  a  >2)(<  62  >  -  <  b  >2) 


Pi  = 


<ab>  -<  ax  b> 
<  a2  >  -  <  a  >2 


Figure  7:  Definition  of  correlation  coefficients  and  an  example  scatter  plot.  In  this 
example  case  for  a  scale-similarity  model  p,  =  0.960  and  p2  =0.519 . 


Results  of  two  a  priori  tests  in  rotating  homogeneous  decaying  turbulence  are  shown  in 
Figure  8.  These  plot  the  correlation  coefficient,  p2 ,  as  a  function  of  the  filter  wave 
number.  The  more  difficult  terms  to  model  are  those  that  have  components  in  the  z 
direction.  Thus,  only  the  stress  term,  r33 ,  and  its  divergence,  dx2jtdxj ,  term  are  shown. 
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The  figures  show  results  for  the  models  presented  here  (SCDSM  (4)-(5)  and  GCDSM 
(15)-(17)).  Results  for  two  other  models  are  also  shown:  DS  is  the  original  dynamic 
structure  model  (8)  and  SSM  is  a  simple  scale  similarity  model  without  rescaling  by  the 
subgrid  kinetic  energy. 

All  models  show  good  results  at  high  wave  numbers  but  are  more  severely  tested  at  low 
wave  numbers  (larger  filter  sizes).  The  results  show  that  the  rescaling  in  the  DS  model 
gives  significantly  better  results  than  the  SSM  model.  As  the  cross  correlation  terms  are 
added  to  the  DS  model  to  give  the  SCDSM  model,  the  correlation  coefficients  improve 
even  more.  The  GCDSM  model  also  gives  good  results  at  larger  filter  sizes  but  goes  to  the 
wrong  limit  as  the  filter  size  decreases. 


k-n/A  k=7i/A 

Figure  8:  Comparison  of  the  second  correlation  coefficient,  p2.  The  filter  cutoff 
wavenumber  is  defined  by  k  =  n/A:  (a)  and  (b)  rotation  rate  fi  =  0.5  [l/s]  and  Rossby 
number  Ro  =  0.91 ,  (c)  and  (d)  rate  Q  =  2.0[l/s]  and  Rossby  number  Ro  =  0.23 .  SSM  is  a 
scale  similarity  model,  DSM  is  the  DS  model  (8),  SCDSM  is  the  dynamic  structure  model 
with  the  cross  terms  (15)-(16)  and  GCDSM  is  a  gradient  model  ( 1 5)-(  1 7). 
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